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Abstract 



We show that a knowledge of diagonal partons at a low scale is sufficient to determine 
the off-diagonal (or skewed) distributions at a higher scale, to a good degree of accu- 
racy. We quantify this observation by presenting results for the evolution of off-diagonal 
distributions from a variety of different inputs. 



1. Introduction 



Precision data are becoming available for hard scattering processes whose description re- 
quires knowledge of off-diagonal (or so-called "skewed") parton distributions |IJ ||. For 
example the diffractive production of vector mesons or high Et jets in high energy electron- 
proton collisions are now experimentally accessible. 

Off-diagonal distributions^ H(x,^,t,fi 2 ) depend on the momentum fractions x\^ = x ± £ 
carried by the emitted and absorbed partons at each scale /i 2 , see the lower part of Fig. 1. They 
also depend on the momentum transfer variable t = (p — p') 2 . However t does not change as 
we evolve the parton distributions up in the scale fi 2 . That is t lies outside the evolution and 
we can study the scale and X{ dependence of H for each fixed value of t. The evolution of the 
distributions H with the scale is described by 



see Fig. 1, where the off-diagonal splitting functions P are known to leading order (LO) and 
the anomalous dimensions are known to next-to-leading order (NLO) accuracy. So there is a 
possibility to evolve H numerically starting from some input distribution. 

As in the diagonal case, we require the form of the starting distributions, which will be 
dependent on non-perturbative QCD. This poses a major problem for off-diagonal evolution 
since we must specify both the x\ and x 2 behaviour of the starting distributions. We could, 
in principle, determine the off-diagonal input directly from data. However the possibility of 
determining the X\ and x 2 behaviour, for a given t value, from data is remote. Alternatively we 
could use some simplified model for the distributions in the non-perturbative region. Clearly 
it would be better to try to relate the off-diagonal distributions to their known diagonal form. 

In the present paper we argue that at a low scale /x 2 = Ql one may input the diagonal forms, 
but smoothed in the ERBL region with \x\ < £ to avoid the singularity at x — > 0. Then after 
a few steps of evolution the result will be very close to the true off-diagonal distribution. In 
other words the original £ dependence is washed out and forgotten during the evolution, with 
the final £ behaviour being generated mainly by the form of the off-diagonal splitting functions. 

To some extent, the idea to use a pure diagonal input is supported by the GRV approach 0. 
With more or less simple, valence-like distributions at very low input scale Qq, the evolution 
generates quite reasonable diagonal parton distributions, whose x-dependence is essentially 
governed by the form of the evolution equation. 

In the next section we give a better argument in favour of diagonal input at a low scale 
for generating reliable off-diagonal distributions by evolution. The argument is based on the 
polynomial properties of the operator product expansion (OPE) for the function H(x,£). 

^^Here we use the off- forward distributions with support — 1 < x < 1 introduced by Ji H, except that 
we take H g = xH^ 1 . In the limit £ — > the distributions reduce to the conventional diagonal distributions: 
H q (x,0) = q(x) for x > 0, H q (x,0) = -q(—x) for x < and H g (x,0) = xg(x). See @, [| || for detailed 
discussion of off-diagonal distributions. 
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In Section 3 we demonstrate numerically that starting at Qq = 0.26 GeV 2 with different 
inputs, which have the same diagonal limit for £ — > 0, we already obtain practically the same 
distribution H(x,£) in the region Q 2 ~ 2 — 10 GeV 2 . Thus one can reasonably specify the 
off-diagonal distributions directly in terms of the (conventional) diagonal partons which are 
known from experiment. 

2. Constraints based on the polynomial condition 

First we recall the known properties^ of the function H(x, £). Just from the Lorentz invari- 
ance and the tensor structure of the operators in the OPE, it follows that the x-moments 

H n (0 = f x n H(x,0dx = £ C kn e k (2) 

are even polynomials in powers of £ of order N < n+1 (i.e. 2k < n+1) @]. The odd powers of 
£ are absent due to the left-right (£ — > — £) symmetry of Fig. 1. Expanding now the off-diagonal 
distribution in the powers of £ 2 

oo 

H(x,0 = J2 ¥ k h2k{x), (3) 

fc=0 

we see that the lower moments of the coefficient functions ti2k{x) with n + 1 < 2k are required 
to be zero 

J x n h 2k (x)dx = 0. (4) 

Therefore the larger the value of k the more the functions h 2 k{x) must change sign and oscillate 
in the interval [—1,1] in order to satisfy the moment conditions. These oscillations will be 
washed out and cancel each other during the evolution. Thus the original (input) £ behaviour 
will die out as Q 2 increases and only the £ dependence generated by evolution will survive. 

Another way to arrive at the same conclusion is to consider the properties of the anomalous 
dimensions 7 n . Here it is better to consider the conformal moments O n |§ 

0* n = J R i n (x 1 ,x 2 ) Hi(x,£) dx, (5) 
with i — q,g, defined on the polynomial bases 

km - ±(" k ) (T + \)**r* ( 6 ) 

- 1 ( i ) ( i ; t ) >?* m 

2 We consider the behaviour at t — 0. Otherwise we may assume the factorization property H(x,^]t, fj, ) = 
G(t)H(x,£; [i 2 ), where G{t) is the proton form factor with G(0) = 1. 
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with xi t 2 = x ± £. The conformal moments have the advantage that they do not mix during 
LO evolution, and are simply renormalized multiplicatively 

O n (Q 2 ) = O n (Ql) (^X\ (8) 

They satisfy the same polynomial properties as the x moments of (0). So only the anomalous 
dimensions 7„ with n + 1 > 2k may participate in the evolution of the coefficient function 
h,2k{x). On the other hand the values of 7 n decrease with increasing n. In fact j n become 
negative for n > 1. Therefore the contributions to H(x, £) coming from the higher terms fi2k{x) 
in @ die out. 

To quantify our observation we perform LO evolution of the off-diagonal distributions from 
a wide variety of different diagonal input distributions. We present a representative selection 
of our results below. To be specific we evolve from different model (M) inputs based on a set 
of diagonal distribution at a low scale Qo, and show, as the scale /i 2 = Q 2 increases, the speed 
with which the models tend to common off-diagonal distributions 

H M (x,H,Q 2 ) -> H(x,t,Q 2 ). (9) 

One possibility is simply to use the diagonal partons themselves as input 

H M (x,Z,Ql) = ho(x) = f{x), (10) 

where f q (x) = q(x) and f g (x) = xg(x). We comment on this below. 

Of course we can add to the above input distributions, functions with support entirely in the 
time-like ERBL region \x\ < £. Evolution never pushes such functions out of the ERBL domain. 
On the other hand as £ — > such contributions become invisible and cannot be generated from 
diagonal parton distributions. However here we discuss functions which are analytic in the 
whole interval — 1 < x < 1. So far the most realistic models for the nucleon (such as fl2fl ) are 
indeed analytic in the whole interval, and do not contain any 0(£ — \x\) contributions. 

For illustration we show results based on two different sets of diagonal partons. In all cases 
we start the off-diagonal evolution from Ql = 0.26 GeV 2 and show results at Q 2 = 4 and 
100 GeV 2 . For set 1 we take the GRV(98) partons and for set 2 the toy model 

xg = 3(1 — a;) 3 , xq = 9x(l — x) 2 , xq = 0, (11) 

which satisfies the momentum and rif = 3 flavour sum rules. Figs. 2 and 3 show the off-diagonal 
distributions obtained from GRV-based inputs at £ = 0.3 and £ = 0.03 respectively, whereas 
Figs. 4 and 5 show the corresponding distributions based on the toy model partons. In each 
plot the dotted curve corresponds to evolution from (unmodified) diagonal partons, as in (|10D . 
We see such an approach gives unphysical input forms in the ERBL region |x| < £, particularly 
for the GRV quark distribution which is singular as x — > 0, that is f q (x) = q(x) ~ x~ a with 
a > 0. It is more realistic to smooth the input in the ERBL region. We show results for two 
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types of smoothing, both of which retain the constraints of the polynomial condition. First the 
continuous curves on Figs. 2-5 are obtained by modifying the diagonal input by requiring the 
conformal moments (||) of H(x,£,Qq) to be independent of £. Second we modify the diagonal 
input by using different forms of h(y) in the double distribution of ref. 0. The dashed curves 
are the result of one interesting choice of h(y). Before we discuss the results shown in Figs. 2-5 
we explain, in turn, the two smoothing procedures. 



In the first case we seek an input which has the diagonal limit H(x,£ = 0; Qjj) = f(x, Q 



a)- 



but whose conformal moments do not depend on £, that is O n {£) = O n (0). This input may be 
calculated by the Shuvaev transform |IUL |TT| 



dx' 
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7T 



i ds (s + e(i-2s))p 

i! y( s )\l l - y(s)x r 



_d_ 

dx' 



(12) 



with p — 1 for gluons and p = for quarks, where 

4s(l 



For small values of £ and x this prescription will in fact generate reliable off-diagonal distri- 
butions from known diagonal partons at any scale ||11|| . For large £ it should produce a more 
physical input than the unmodified diagonal form (|T0|) , and will much more quickly evolve into 
the true off-diagonal distributions. The results are shown by the continuous curves in Figs. 2-5; 
the smoothing of the distributions in the ERBL-region is manifest. 

An alternative way to smooth out the input distributions is to use the double distributions 
of ref. || 

H(x,0 = t F(x-Cy,y)dy (13) 

J a 

with 



F{x,y) = f{x) h(y) 



l-x 



h(y') dy 



where h(y) is an even function of y and f(x) is a diagonal distribution. Thus in the limit £ 
we have H(x,$, = 0) = f(x). The integration limits in ( |T3"D are 



(14) 
= 



[a, b) 
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for x > £ 

for — ^ < x < £ 



(15) 



This construction satisfies the polynomial condition. Indeed, if we write relation (|I3]) in the 
equivalent form 

H(x,0 = I dxdyF(x, y )8(x-(x + Zy)) (16) 



4 



where the integration region R is given by the square in the (x, y) plane defined by the vertices: 
(0, ±1) and (±1,0), then condition (0) is explicitly satisfied. 

We are free to choose the form of h(y) and generate a whole range of input distributions. 
The choice 

%) = (l-y 2 Y + \ (17) 

with p — 0, 1 for quark (gluons), generates an input form which turns out remarkably similar to 
that obtained by the Shuvaev transform (|T2|) . We tried several other forms for h(y). A relevant 
choice is 

h(y) = sm(ny 2 ) (18) 
for both quarks and gluons, since it generates an oscillatory input behaviour in the ERBL 



region which is qualitatively similar to the quark distribution calculated in ref. p2[ . However 
after evolving to Q 2 = 4 GeV 2 we see that the oscillations are much weaker, and that by 
Q 2 = 100 GeV 2 they have already disappeared. 

The similarity of the input based on dl2|) and fll7l) at scale Ql = 0.26 GeV 2 is fortuitous. 
The Shuvaev transform (|i~2"D has the advantage that it is based on conformal moments which 
renormalize multiplicatively, (^), with the same anomalous dimensions 7 n as in the diagonal 
case. So the transform is stable to a change of scale, whereas this is not the case for fllTD . 
Therefore prescriptions based on ([TJ) and (pj]) will give different off-diagonal distributions if 
applied at a higher scale. 

From the representative results shown in Figs. 2-5, and from those obtained by other trial 
inputs based on diagonal distributions, we can draw the following conclusions. 

(i) If we are given a set of diagonal partons f(x, Qq) then we can predict the off-diagonal 
distributions H(x, £, Q 2 ) at a higher scale with a high degree of certainty. Any differences, 
between the prescriptions used, rapidly disappear with increasing Q 2 . 



(ii) The prescriptions are most precise at small values of £, as expected from ref. [11. 



(iii) There is freedom in the form of the input distributions in the ERBL region, |x| < £, but 
even here the various choices evolve to a common limit reasonably rapidly with increasing 
Q 2 , as shown, for example, by comparing the continuous and dashed curves in Figs. 2-5. 
(The dotted curves, corresponding to unmodified input partons, are shown simply for 
comparison.) 

(iv) For gluons in the whole domain, and for quarks in the DGLAP region \x\ > £, the 
distributions evolve particularly quickly to a common limit. 

These conclusions are encouraging for phenomenology. For instance they will enable data for 
processes described by off-diagonal distributions to provide additional independent constraints 
in the global (diagonal) parton analyses. For example, one type of process for which data are 
becoming available is high energy diffractive vector meson electroproduction. The description 
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is dependent on the off-diagonal gluon distribution at small x, £, which is especially well known 
in terms of the diagonal gluon distribution, see Figs. 3 and 5. 
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Figure 1: A schematic diagram showing the variables for the off-diagonal parton distribution 
H(x,£), where — x ± £. The timelike region with \x\ < £ is frequently called the ERBL 
domain, since in the limit £ — ■> f "pure" ERBL evolution applies. Similarly \x\ > £ is known as 
the DGLAP domain since we have "pure" DGLAP evolution for £ — > 0. 
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GRV 5=0.3 





Figure 2: The evolution to Q 2 = 4 and 100 GeV 2 of the quark non-singlet H ns , the quark 
singlet H s and the gluon distribution H g (= xH^ 1 ) starting from diagonal GRV-based input at 
Ql = 0.26 GeV 2 for £ = 0.3. The dotted curves correspond to evolution from unmodified input 
partons as in (0), the continuous curves to using the "Shuvaev prescription" of (|12D and the 
dashed curves to "double distribution" input with (|18|). 
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GRV ^=0.03 




Figure 3: As in Fig. 2 but for £ = 0.03. 
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Toy model £,=0.3 




Toy model ^=0.03 




Figure 5: As in Fig. 4 but for £ = 0.03. 
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